Nonlinear spin-up of a thermally stratified fluid in 
cylindrical geometries 



J. Rafael Pacheco^'^, Sergey A. Smirnov^, and Roberto Verzicco'^ 

^Arizona State University, USA 
^The University of Notre Dame, USA 
^Texas Tech. University, USA 
^Universita' di Roma "Tor Vergata" , Italy 

October 26, 2009 



Abstract 

This is an entry for the Gallery of Fluid Motion of the 62nd An- 
nual Meeting of the APS-DFD (fluid dynamics videos). This video 
shows the three-dimensional time-dependent incremental spin-up of a 
thermally stratified fiuid in a cylinder and in an annulus. The rigid bot- 
tom/side wall(s) are non-slip, and the upper surface is stress-free. All 
the surfaces are thermally insulated. The working fluid is water charac- 
terized by the kinematic viscosity v and thermal diffusivity n. Initially, 
the fluid temperature varies linearly with height and is characterized 
by a constant buoyancy frequency A'^, which is proportional to the den- 
sity gradient. The system undergoes an abrupt change in the rotation 
rate from its initial value Qj, when the fluid is in a solid-body rotation 
state, to the final value fij. Our study reveals a feasibility for transi- 
tion from an axisymmetric initial circulation to non-axisymmetric fiow 
patterns at late spin-up times. 

1 Introduction 

This fluid dynamics video shows the baroclinic instability of the flow evolv- 
ing in time due to an abrupt change of the rotation rate of the cylin- 
der/annulus from diff^erent perspectives in space. Our simulations revealed 
that azimuthal asymmetry is manifested in the form of cyclonic and anti- 
cyclonic columnar eddies, which develop at the temperature front formed 



by the highly distorted isotherms near the cyhnder sidewalls. Strong defor- 
mation of the initial temperature field is caused by the Ekman transport 
near the bottom. The front steepens until it reaches a quasi-equilibrium 
state. The eddies grow in size and march along the outer wall until they 
occupy a large portion of the tank. The characteristic Ekman number E 
is small and the Ekman bottom boundary layer is in the transitional state. 
The Reynolds number based on the Ekman layer depth is Res = 90 and the 
Rossby number e = {Qj — Qi)/Qf (=0.727), so that the flow regime under 
investigation is highly nonlinear. The details of the numerical method and 
validation tests are described in [3] and references therein. It is important 
to note that the formation of these large-scale eddies has been observed ex- 
periments of stratified fluids by salt [H El [2]. The variables shown in the 
video are the temperature contours and isosurfaces of the Q-criterion [T] 
colored by temperature. The first section of the video presents the three- 
dimensional numerical simulations for the cylinder, and the second section 
of the video corresponds to the spin-up in the annulus. 
The video animation file can be found at spinup.mpg] 
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